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ABSTRACT 

We examine the constraints imposed on the total kinetic power (Lj) and the age 
(t ag e) of relativistic jets in FR II radio sources in a new way. We solve the dynamical 
expansion of its cocoon embedded in the intra-cluster medium (ICM) and obtain the 
analytic solution of its physical quantities in terms of Lj and i age . The estimate of 
Lj and t agc is done by the comparison of the model and the observed shape of the 
cocoon. The analysis is focused on Cygnus A and we find that (i) the source age is 
3 Myr < t agc < 30 Myr and (ii) the total kinetic power of the jet is estimated as 
0.2 x 10 46 erg s _1 < Lj < 1 x 10 48 erg s^ 1 which is larger than 1% of the Eddington 
luminosity of Cygnus A. 

Key words: plasmas — radiation mechanisms: non-thermal — galaxies: active — 
galaxies: individual: Cygnus A 



1 INTRODUCTION 

After the detections of inverse Compton component in X- 
ray from FR II radio sources, our knowledge of energetics, 
especially on the kinetic power of non-thermal electrons, is 
progressed in recent years (e.g., Leahy and Gizani 2001; Har- 
ris and Krawczynski 2002; Hardcastle at al. 2002; Kataoka 
et al. 2003; Croston et al. 2005). When we explore further 
physical conditions in the jets, we encounter a big difficulty 
to constrain on the fraction such as thermal electrons and/or 
protons co-existing with non-thermal electrons because it is 
hard to observe these component (e.g., Celotti et al. 1998; 
Sikora and Madejski 2002). This problem prevented us from 
estimating the total mass and energy flux ejected from a cen- 
tral engine. To conquer this, a simple procedure is proposed 
in the study of jets in FR II sources in Kino and Takahara 
(2004) (hereafter KT04). In the pressure and mass density 
of jet, the contributions from the invisible components are 
also involved. Hence the rest mass density and energy den- 
sity estimated from the shock dynamics definitely prove the 
quantities of total plasma. Also for radio bubbles in cluster 
cores, the similar dynamical approach has been adopted to 
constrain on their physical state (Fabian et al. 2002; Dunn 
and Fabian 2004). 

In this paper, we will explore the cocoon dynamics in 
FR II radio sources as a robust tool to know the total ki- 
netic power Lj and source age t ag o. In §2, we discuss our 
analytical solution of cocoon expansion which includes the 
effect of radial dependence of the rest mass density of the 
surrounding Intra Cluster Medium (ICM). In §3, by com- 



paring the analytic solution and observation of the cocoon, 
we explore the kinetic power of Cygnus A, which is one of 
the best-studied FR II source. Conclusions are given in §4. 



2 DYNAMICS OF COCOON 

We deal with the cocoon dynamics in FR II radio source. 
The adopted basic equations are almost the same as those 
in Begelman and Cioffi (1989) (hereafter BC89). The main 
differences between BC89 and the present work are (i) we 
explicitly solve the physical quantities as a function of Lj, 
and f ago , (ii) we take account of the effect of radial profile of 
the mass density of ICM, and, as a result, (iii) the growth 
of the cocoon head Ah is consistently solved from the basic 
equations. 

2.1 Basic Assumptions 

Our main assumptions are as follows; (1) We limit our at- 
tention a jet relativistic speed (vj = c), (2) The jet supplies 
a constant Lj in time, (3) We focus on the over-pressured 
cocoon phase, and (4) We assume that the magnetic fields 
are passive and ignore their dynamical effect. For (1), al- 
though it is still under debate, some jets are suggested to be 
relativistic (e.g., Tavecchio et al. 2000; Celotti, Ghisellini, & 
Chiaberge 2001). Since little is known about the evolution 
of Lj, the assumption (2) is adopted as a first-step working 
hypothesis. The assumption of (3) is automatically guaran- 
teed by the sideways expansion of the cocoon (e.g., Cioffi 
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and Blondin 1992). The assumption of (4) is based on the 
results that a multi-frequency analysis of radio galaxies show 
that the energy density of magnetic field tend to be smaller 
than that of non-thermal electrons (e.g., Leahy and Giani 
2001; Isobe et al. 2002). 

2.2 Basic Equations 

Model parameters are Lj and i agc . Unknown physical quan- 
tities are Vh, v c , Pc, and Ah (or A c ) which are the advance 
velocity of the cocoon head, the velocity of cocoon sideways 
expansion, the pressure of the cocoon, the cross sectional 
area of the head part of the cocoon (or the cross sectional 
area of cocoon body), respectively (see Fig. 1). Equation 
of motion along the jet axis, and sideways expansion, and 
energy conservation in the cocoon are given by 

U 



Pa.(rh)v h (t)A h (t), 



P c (t) = p a (r c ) v c (t) 2 



Pc(t)V c (t) 



2L s t, 
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where r h (t) = J^v h (t')dt', r c (t) = f^v c (t')dt', V c (t) = 
2 J A c (t')vh{t')dt' , t m jn, and 7 C are the length from the 
center of the galaxy to the head of the cocoon, the radius of 
the cocoon body, the volume of the cocoon, the start time 
of source evolution and specific heat ratio of the plasma 
inside the cocoon, respectively. The declining mass density 
of ICM p a is given by p a (r) = p a (r/ro) _c * where ro and 
p a are reference position and the ICM mass density at ro, 
respectively. We set ro as rh(£ agc ) in throughout this paper. 
Most of the kinetic energy is deposited in the cocoon, which 
is initially suggested by Scheuer (1974) and recent studies 
of hot spots also shows the radiative efficiency is very small 
(e.g., KT04). Following to Cioffi and Blondin (1992), we add 
the factor of in Eq. to express the amount of the 

deposited internal energy. In other words, this corresponds 
to the neglect of the PdV work because of its smallness. 

The numbers of quantities are 4, while those of basic 
Eqs. are 3. Hence, we set A c (t) oc t x where X as a free 
parameter. We can constrain on the value of X from obser- 
vations. A specific case is shown in §3. Hence, we obtain Vh, 
v c , P c , and Ah by using a free parameter X. 

As a subsidiary equation, the area of the cocoon body 

is given by A c (t) = tt (^f* v c (t')dt'^J . It is clear that the 

change of the unknown from A c to Ah does not change the 
result. However note that when we choose Ah as an unknown 
instead, we cannot obtain the solution for a = 2. 

2.3 Analytic Solution 

We assume that physical quantities have a power law de- 
pendence in time and the coefficient of each physical quan- 
tity is barred quantity. Each quantity has the form of 
A = A (t/t age) where Y is an arbitrary index. The time 
evolution of v c is 



v c (t) 



t 
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for a given A c . With this, the analytic form of cocoon quan- 
tities in decreasing ICM density is obtained as follows; 
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Ah{Lj,t) = 



Li 
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where C vh = (•% - 1)[3 - (1 - 0.5«)X](0.5X)- a , C vc = 
0.5X/tv 1/2 , C pc = (0.5A) Q , andC ah = [X(-2 + 0.5a) +3]" Q , 
respectively, vo = rh/t ag e corresponds to the head speed for 
constant velocity in time. 

Here we use the conditions of 0.5X > and X(— 2 + 
0.5a) + 3 > which make the contribution at t mln in the 
integrations of rh and r c small enough. The case we focus on 
in §3 is that X = 12/7 and a — 1.5, which clearly satisfies 
these conditions. 

First, let us consider the evolution of cocoon. The 
growth of both Ah and A c must be positive. As for the co- 
coon expansion speeds, three different behaviors are theoret- 
ically predicted such as (I) accelerated-head (X( — 2 + 0.5a) + 
2 > 0), (II) constant-head (X(-2 + 0.5a) + 2 = 0), and (III) 
decelerated-head (X(-2 + 0.5a) + 2 < 0). The case of (I), 
(II), and (III) correspond to X < 1, X = 1, and X > 1 for 
a — 0, while in the case of a = 2 (I), (II), and (III) corre- 
spond to X < 2, X — 2, and X > 2, respectively. Related 
to this, it is useful to express the aspect ratio of the cocoon 
— = 1Z as a function of time. This is written as 
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It is worth to note that the solution describes not only the 
self-similar evolution (e.g., Kaiser and Alexander 1997; Bick- 
nell, Dopita and O'Dea 1997) but also the non self-similar 
one. Although it is not concerned with in this paper, the 
evolution of lZ(t) may probe the evolution of radio galax- 
ies such as "compact symmetric objects (CSO)" which are 
thought to be the progenitors of the FR II source (e.g., Fanti 
et al. 1996; Readhead et al. 1996a). Note that, observation- 
ally, the large deviation from 1Z ~ 0(1) does not seem to 
be unnatural for actual radio sources (e.g., Readhead et al. 
1996b). 

Next, let us consider the a and X dependences on the 
quantities based on Eqs. (1), (2) and (3). Here we fix the 
physical quantities at t — i a gc- With regard to the a de- 
pendence in fixed X, larger a leads to larger p a , slower Vh, 
same v c , larger P c , during t < t ag e. We can understand it 
as follows. Larger a leads to stronger the deceleration effect 
on the head speed Vh due to larger p a . Larger a also predict 
larger P c in order to keep the same velocity of sideways ex- 
pansion v c - In Fig. H we show the effect of varying a when 
fixing X = 1. It shows that larger a requires larger Lj in 
order to plow the larger amount of ICM. 

Next we consider the X dependence. In fixed a, larger 
X lead to faster Vh, slower v c , smaller P c , and smaller Ah 
during t < t agc . We can explain it as follows. Lager X leads 
to slower sideways expansion v c and smaller A c by definition. 
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From the equation of motion to the sideways expansion, it 
is clear that smaller P c is required. To satisfy the energy 
equation at the same time, the faster velocity of Vh is needed. 
This is realized by smaller area head of the cocoon Ah. 

Observationally, little is argued about the emission from 
the cocoon itself. Recently, Readhead et al. (1996b) report 
the cocoon emission from the CSO 2352+495. It is visible in 
the 610 MHz image and the aspect-ratio is about H(t m in) ~ 
1/2. Because of the lack of our knowledge of the cocoon 
emissions, here we examine the case of 7?.(i m in) =1, 1/2, 
and 1/4 for the moment. In Fig. [3] we show the effect of 
varying X with fixing a = 1.5 and TZ(t me ) = 1/2. The case 
of X = 1.62, 12/7, and 1.80 correspond to the K(t mln ) =1, 
1/2, and 1/4 respectively. We can verify the shift of Lj to 
smaller range as the X increases. 

Finally, we compare our solution with the previous 
works. The solution for flat ICM density by BC89 model 
corresponds to the case of X = 1 and a — and the quan- 
tities are written as v c oc t^ 1 ^ 2 , Vh oc const., P c oc t" 1 , 
Ah oc const. The absolute values of each physical quantities 
are also agree with those in BC89 when omitting — 1) 

and replace C vc = I/tt 0,5 , C v h = 1, C pc = 1, and C a h = 1. It 
is included in Fig. [5] the case of BC89 with the parameters 
of a = 0, X = 1, and Ah — 30 kpc 2 . As explained before, 
the larger (smaller) a predict smaller (larger) Lj because of 
the dilute (dense) ICM. The comparison with the analytic 
model and numerical studies, is also important issue. By 
Cioffi and Blondin (1992), the issues of "head cross section 
growth" and "decreasing head velocity" is assessed by hy- 
drodynamic simulation and their result shows Ah oc t 0A . As 
a good example, our solution of X ~ 1.1 describes the stud- 
ies Nath (1995) corresponds to case of flat solution above 
Vh oc t~ - 2 , v c oc i~ ' 45 , P c oc i~ ' 9 , and A h oc t 0A . Concern- 
ing A c , since most of numerical studies have mainly focused 
on the propagation of cylindrical geometry jets (e.g., Marti 
et al. 1997; Clarke, Harris and Carilli 1997), it remains as 
an important future work on the time evolution of A c . 



/ f \ -11/7 

Pc(t) = 4.78 x lO^'Vo 2 ( — d y ne cm ~ 2 > ( n ) 

A b (t) = 66.4 L\%t\(? (±"j kpc 2 , (12) 

where f3 h = v h /c, C = v c /c, t 2 o = t agc /20Myr, L jAe = 
Lj/10 46 erg s , with the resultant value of 7Z(t agc ) = 0.815 
for Cygnus A. Note that TZ oc /3h/Pc has a strong dependence 
on the age such as oc t 3 go , the allowed source age is fairly 
restricted. 

Next, to clarify the allowed range of t ago and Lj, we 
impose the following conditions; the condition which should 
be satisfied is that (I) 7£(t ag c) ~ 0.5 — 0.7, from the Chandra 
image (Wilson et al. 2000), (II) cocoon pressure is over- 
pressured P c > P& = 8x 10 _11 dyn cm -2 (Arnaud et al. 1984; 
Smith et al. 2002) , (III) the area size of Ah lies in the range of 
30 kpc 2 < A h < 150 kpc 2 . The minimum value corresponds 
to the one adopted in BC89. From the radio image of Perley, 
Dreher and Cowan (1984), we employ the maximum value as 
Ah = 150 kpc 2 , which corresponds to the cross sectional area 
of the radio lobe at the location of the hot spot. In Fig. 01 we 
show the resultant source age and total kinetic power of the 
jet. The region of the source age larger than ~ 30Myr is not 
allowed by the condition (II). Larger (smaller) Ah predict 
larger (smaller) Lj and younger i agc . Obtained values are 
2 x 10 45 erg s" 1 < Lj < 1 x 10 48 erg s" 1 and 3 Myr < i agc < 
30 Myr. 

Let us consider how the uncertainties of Lj and i age 
are determined. From Eqs. Q and JHJ, Ah oc Lji agc and 
TZ oc LT~ 1 ^ ct ~ 4 '^t^gJ ( ' a are obtained. Since Ah and TZ have 
uncertainties with the factors of 5 and 1.4 respectively, the 
allowed region is mainly controlled by Ah- In the present 
work, in view of the comparison of BC89, we took the min- 
imum value as Ah = 30kpc 2 . However, from the hydrody- 
namical point of view, it seems natural to suppose that Ah 
to be the close value to the cross sectional area of the radio 
lobe at the distance of head, i.e., Ah = 150 kpc 2 . 



3 TOTAL KINETIC POWER AND SOURCE 
AGE OF CYGNUS A 

Here we explore Lj and £ ago by matching the observed cross- 
sectional areas and lengths of a cocoon (i.e., rh, r c , Ah, and 
A c ). In this paper, we focus on the archetypal radio galaxy 
Cygnus A (Carilli and Barthel 1996; Carilli and Harris 1996 
for reviews). 

First, we estimate the typical values of physical quanti- 
ties. Concerning the mass density profile of ICM, we adopt 
q = 1.5 based on Reynolds and Fabian (1996) and Smith et 
al. (2002). We examine the case of X = 12/7 as a fiducial one 
which predict the constant TZ in time. Other observed quan- 
tities p a = 0.5 x 10 _2 m p g cm -3 , and rh = 60 kpc, based 
on Carilli et al. (1998), Arnaud et al.(1984), and Smith et 
al. (2002). Here we assume 7 C = 4/3. Using these quantities, 
we obtain 

f3 h (t) = 8.36 x 10" 3 Lj, 46 4) (^) 1A , (9) 
&(*) = 6.84 X 10" 3 t^ 1 f-M 1A , (10) 



3.1 Comparison with previous works 

The resultant age well agree with the independent result of 
synchrotron age model by Carilli et al. (1991), which claims 
that 6 Myr < t agc < 30 Myr. The velocity of the hot spot 
Phs ~ 0.06 corresponds to the source age of 6 Myr, while 
Phs ~ 0.01 corresponds to the source age of 30 Myr. 

As a complementary result, Lj estimated in KT04 in 
the range of 6 Myr < t ago < 30 Myr based on the result 
Carilli et al. (1991) is shown in Fig. ||| (the solid line). KT04 
estimate the total kinetic power of the relativistic jet as 

I- = l '><•'' 

= Ai c ( 7^" ) Pa(r 60 ) oc Aj 
yt ag c J 

in the strong relativistic shock limit, where Aj = ttP 2 s; Fj, 
/3jC, pj, and Phs = 2 kpc (Wilson et al. 2000) are the cross- 
sectional area of the jet the Lorentz factor, the velocity, and 
mass density of the jet, and the hot spot radius, respectively. 
It should be stressed that the cocoon model can predict both 
Lj and t age at the same time, while the synchrotron aging 
model (Carilli et al. 1991) and the ID jet model (KT04) only 
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determine i age or Lj . Although these three models are inde- 
pendent, they show a reasonable agreement on the values of 
Lj and t ago with each other at least in order- of- magnitude 
in the case of Cygnus A. At the same time, it is worth to 
discuss a factor of deviation of the estimate Lj by the co- 
coon model and ID jet model. In the same way as KT04, 
ID analysis between the jet and ICM take the plane parallel 
assumption (e.g., Begelman, Blandford & Rees 1984). The 
approximation means that Aj equals to A^. However, A^ 
is supposed to be larger than Aj. Hence, the plane parallel 
approximation would cause the underestimate of the Lj in 
KT04. 

Rawlings & Saunders (1991) (hereafter RS91) is the pi- 
oneering work on the correlation between the Lj and the 
luminosity of the narrow line regions, which lies in close to 
the central engine. Then the comparison with the present 
work and RS91 is intriguing issue. For the sources where 
synchrotron spectral aging are available, it is simply by 



where Q, rj, and E ecl are the kinetic power of jet, a parameter 
expressing the fraction of the work done on the ICM, and 
the equipartition energy with the electrons and the magnetic 
field field make an equal contribution to the total energy 
density (e.g., Miley (1980)), respectively. We focus on the 
sample sources where i age is independently obtained by syn- 
chrotron aging method. One difference between the present 
work and RS91 is that we solved the equations of motion 
and energy equation (3 Eqs. in total), while the RS91 only 
employ energy equation. Because of this, we can eliminate 
the free parameter r). More important and essential differ- 
ence is that Q in RS91 is not a total kinetic power but it 
is just a equipartition power (even though they insist that 
it as a total power). We emphasize the advantage of our 
work of dealing with the total kinetic power whilst RS91 
only handles the equipartition power of extragalactic jets. 

Kaiser (2000) (hereafter K00) independently addressed 
this quantity by studying some bright radio sources involv- 
ing Cygnus A. Compared with RS91, the common advantage 
is that both K00 and the present work can estimate the total 
kinetic power by directly dealing with the hydrodynamics. 
The way of comparing the observation with the model by 
K00 is different from ours. K00 matched the surface bright- 
ness distribution of the cocoon along the jet axis. The main 
difference between the present work and K00 is the derived 
cocoon pressure during the matching of the observations and 
the models. In the example of Cygnus A, compared with our 
estimate of P c ~ 5 x 10~ 10 erg cm -3 , they tend to estimate 
smaller P c of order of P c ~ 10~ 12 - 10~ n erg cm" 3 (Tables 
2, 3 and 4 in K00). This mainly cause the deviation of the de- 
rived total kinetic power of Cygnus A as Lj ~ (a few) x 10 46 
erg s~ , whilst K00 derive the total kinetic power with order 
of Lj ~ (a few) x 10 45 erg s _1 . 

3.2 Efficiency of accretion power to kinetic power 

In Tadhunter et al. (2003), the mass of the super-massive 
black hole (SMBH) of Cygnus A is reported as 2.5 x 10 9 M© 
which leads to the Eddington luminosity as L E dd = 3 x 10 47 
erg s _1 . From this it follows that Lj/L B dd ~ 0.01 — 1. On the 
basis of the observational evidence for the accretion-flow ori- 



gin for the jet in AGNs (e.g., Marscher et al. 2002), it is clear 
that the value of Lj /LEdd directly shows the required mini- 
mum rate of the mass accretion onto the SMBH normalized 
by the corresponding Eddington mass accretion rate. 

Merloni, Heinz and Di Matteo 2003 (see also Mac- 
carone, Gallo and Fender 2003) examine the disk-jet con- 
nection by studying the correlation between the radio (Lr) 
and the X-ray (Lx) luminosity and the black hole mass. 
With the large samples of the stellar mass black holes and 
SMBHs, they claim these sources have the "fundamental 
plane" in the 3 dimensional space of these quantities. On 
the way, the quantity of Lx /LEdd is brought up to probe 
the activity of the central engine. For Cygnus A, Young et 
al. (2002) reports L x ^ 3.7 x 10 44 erg s _1 which leads to 
Lx/LEdd ~ 10 . Whereas we recognize the usability of 
the quantity of Lx/L E dd to probe the engine activity, we 
emphasize that the Lx largely depends on the accretion 
flow model and the radiation efficiency at X-ray band. On 
the other hand, the quantity of Lj/Lsdd addressed in the 
present work does not have the model dependence and di- 
rectly shows us the engine activity. 

Maraschi & Tavecchio (2003) also addressed this quan- 
tity by collecting large samples of blazars. They estimated 
the jet power from the blazar's non-thermal spectrum en- 
ergy distribution, which are well understood as synchrotron 
plus inverse Compton (e.g., Maraschi et al. 1992; Sikora et 
al. 1994; Blandford & Levinson 1995; Kino et al. 2002). The 
upper limit of the accretion power is inferred from the lumi- 
nosity of observed broad emission lines. They conclude that 
for flat-spectrum radio quasars (FSRQs) the the total power 
of the jet is of the same order as the accretion power. Note 
that their estimate of the accretion power from the emission 
line is putative way and the estimated total power does not 
include the contribution of protons associated with thermal 
electrons which may cause some underestimate of the total 
power of the jet. 

We again emphasize that as we see above, the value of 
Lj/Lsdd for Cygnus A in the present work is obtained with 
fewer assumptions. Moreover it make us a fairly robust probe 
of the central engine activity of AGN jets. The application 
of our method to a large sample of other AGN jets surely 
bring about new and important knowledges. This is actually 
our ongoing project. 



4 CONCLUSION 

Our main conclusions in the present work are as follows; 

(i) A new method to estimate the total kinetic power 
of the jet Lj and source age t agc in powerful FR II radio 
sources is proposed. For that, the study of cocoon dynamics 
by BC89 is revisited and physical quantities associated with 
the cocoon are analytically solved as functions of Lj and f agc . 
The comparison of the analytic solution with the observed 
cocoon shape lead to the Lj and t age in general. 

(ii) The analysis is focused on Cygnus A, with the con- 
ditions of 0.5 < TZ < 0.7 and 30 kpc 2 < A h < 150 kpc 2 . 
The estimated age 3 Myr < i agc < 30 Myr shows a good 
agreement of the independently estimate age by synchrotron 
aging model by Carilli et al. (1991). The estimated total ki- 
netic power 0.2 x 10 46 erg s _1 < Lj < 1 x 10 48 erg s _1 has 
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also a reasonable agreement with the independent ID jet 
model of KT04 with the aid of 6 Myr < i agc < 30 Myr. 

(iii) For Cygnus A, we find that the total kinetic power 
lies in the range of Lj/Lsaa ~ 0.01 — 1, while the X-ray lu- 
minosity Lx — 3.7 x 10 44 erg s _1 (Young et al. 2002) satisfies 
Lx/LEdd ~ 10~ 3 . The result of Lj/Z/Edd ~ 0.01 — 1 indicate 
the lower limit of the mass accretion rate, which gives the 
crucial hint for resolving the jet formation problem. 
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Figure 2. The allowed regions of Lj and i age with the differ- 
ent values of a with X 1. The region is determined by the 
constraints of 0.5 < TZ < 1 and 30kpc 2 < A h < 150kpc 2 . We 
examine a =0 and 2 here. The case of a = 0, and = 30kpc 2 
by Begelman and CiofB (1989) is involved in the right-lower part 
of the filled region. Larger a requires significantly larger Lj for 
plowing the larger amount of ICM. 
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Figure 4. The allowed region of Lj and f agG of Cygnus A (filled 
in black) . The under-pressured region P c < P a is excluded by def- 
inition. The case of 0.5 < TZ < 0.7 and 30kpc 2 < A h < 150kpc 2 is 
examined. As a reference, Eddington luminosity and the total ki- 
netic power of jet estimated in KT04 are shown in the thick-solid 
and solid lines, respectively. 
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Figure 3. The allowed regions of Lj and t agc with the diffferent 
values of X with a = 1.5. Likewise Fig. |2] we adopt 0.5 < TZ < 1 
and 30kpc 2 < A h < 150kpc 2 . The case of X =1.62, 12/7 and 
1.80 are shown here. Lager X requires smaller Lj corresponding 
to the slower velocity of the sideways expand. 



